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We show compatibility with all known experimental constraints of t — b — r Yukawa coupling unification 
in supersymmetric SU(4) C x SU(2)l x SU(2)r which has non-universal gaugino masses and the MSSM 
parameter ji < 0. In particular, the relic neutralino abundance satisfies the WMAP bounds and A(g — 2) M is 
in good agreement with the observations. We identify benchmark points for the sparticle spectra which can be 
tested at the LHC, including those associated with gluino and stau coannihilation channels, mixed bino-Higgsino 
state and the A-funnel region. We also briefly discuss prospects for testing Yukawa unification with the ongoing 
and planned direct detection experiments. 



, I. INTRODUCTION 

r-H 

Supersymmetric 50 (10), in contrast to its non- supersymmetric version, yields third family (t — b — r) Yukawa unification via 
the unique renormalizable Yukawa coupling 16-16-10, where the 10-plet is assumed to contain the two minimal supersymmteric 
standard model (MSSM) Higgs doublets H u and The 16-plet contains the 15 chiral fermions per family of the standard 
model (SM) as well as right handed neutrino. The implications of this unification have been extensively explored over the years 
[1, 2]. More recently, it has been argued in [3, 4] that 50(10) Yukawa unification predicts relatively light (< TeV) gluinos, 
which can be readily tested [5] at the Large Hadron Collider (LHC). The squarks and sleptons turn out to have masses in the 
multi-TeV range. Moreover, it is argued in [3, 4] that the lightest neutralino is not a viable cold dark matter candidate, at least in 
i the simplest models of 50(10) Yukawa unification. 

Spurred by these developments we have investigated t — b — r Yukawa unification [4, 6] in the framework of supersymmetric 
Qh SU(4) c x SU(2)l x SU(2)r [7] (4-2-2, for short), with a positive supersymmetric bilinear Higgs parameter (/i > 0). The 
4-2-2 structure allows us to consider non-universal gaugino masses while retaining Yukawa unification. In particular, assuming 
left-right symmetry [7, 8], or more precisely C-parity [9], we have only one additional free parameter in the soft supersymmetry 
breaking (SSB) sector compared to the 5O(10) model. In particular, this allow us to distinguish M2 from M3, where M2 
(Ms) denotes the asymptotic SU (2) l (SU (3) c ) gaugino mass. An important conclusion reached in [4, 6] is that with gaugino 
non-universality, Yukawa unification in 4-2-2 is compatible with neutralino dark matter, with gluino co- annihilation [4, 6, 10] 
^ playing an important role. 

The main purpose of this paper to extend the 4-2-2 discussion to the case of \i < 0, where /i denotes the supersymmetric Higgs 
s O mass parameter. Most authors normally do not consider fi < because it can create serious disagreement with the measured 
value of the muon anomalous magnetic moment (g — 2) M [11]. The new contribution to (g — 2) M from supersymmetric particles 
is proportional to /iM 2 tan f3/m 4 [12], where fh is the heaviest sparticle mass in the loop. In 5O(10) Yukawa unification, the 
scalar masses are very heavy (multi-TeV), so that the supersymmetry (SUSY) contribution to (g — 2) M effectively decouples and 
one is left with the standard model result [3, 4]. On the contrary, in 4-2-2, with both fi < and M 2 < 0, the SUSY contributions 
to (g — 2) M turn out to be important, as we will see, and thereby yield significantly improved agreement with experimental data. 
The outline for the rest of the paper is as follows. In Section II we briefly describe the model and the boundary conditions 
^ for SSB parameters which we employ for our scan. In Section III we summarize the scanning procedure and the experimental 
constraints that we have employed. In Section IV we discuss how /i < leads to better Yukawa unification than fi > 0. 
In Section V we present the results from our scan and highlight some of the predictions of the 4-2-2 model. The correlation 
between the spin-independent and spin-dependent direct detection of dark matter and Yukawa unification condition is presented 
in Section VI where we also display some benchmark points. Our conclusions are summarized in Section VII. 



II. THE 4-2-2 MODEL 



In 4-2-2 the 16-plet of 5O(10) matter fields consists of xp (4, 2, 1) and %jj c (4, 1, 2). The third family Yukawa coupling ipcipH, 
where H( 1,2,2) denotes the bi-doublet (1,2,2), yields the following relation valid at Mqut, 

Y t =Y h = Y T = Y VT . (1) 
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Supplementing 4-2-2 with a discrete left-right (LR) symmetry [7, 8] (more precisely C-parity) [9] reduces the number of 
independent gauge couplings in 4-2-2 from three to two. This is because C-parity imposes the gauge coupling unification 
condition (g^ = gR) at Mqut- We will assume that due to C-parity the SSB mass terms, induced at Mqut through gravity 
mediated supersymmetry breaking [13] are equal in magnitude for the squarks and sleptons of the three families. The tree 
level asymptotic MSSM gaugino SSB masses, on the other hand, can be non-universal from the following consideration. From 
C-parity, we can expect that the gaugino masses at Mqut associated with SU (2) l and SU (2)r are the same (M 2 = M 2 R — 
M2). However, the asymptotic SU (4) c and consequently SU (3) c gaugino SSB masses can be different. With the hypercharge 
generator in 4-2-2 given by Y = y/2/5 (B — L) + \/^J^> Isr, where B — L and I%r are the diagonal generators of SU (4) c and 
SU (2)r, we have the following asymptotic relation between the three MSSM gaugino SSB masses: 

M 1 = ^M 2 + ?M 3 . (2) 

The supersymmetric 4-2-2 model with C-parity thus has two independent parameters (M 2 and M 3 ) in the gaugino sector. In 
order to implement Yukawa unification it turns out that the SSB Higgs mass terms must be non-universal at Mqut- Namely, 
m H u < 171 H d at Mgut, where mriH u ( m H d ) is the up (down) type SSB Higgs mass term. The fundamental parameters of the 
4-2-2 model that we consider are as follows: 



m , m Hu , m Hd ,M 2 , M 3 , A , tan /3, sign(/i). (3) 

Here mo is the universal SSB mass for MSSM sfermions, A is the universal SSB trilinear scalar interaction (with the corre- 
sponding Yukawa coupling factored out), tan f3 is the ratio of the vacuum expectation values (VEVs) of the two MSSM Higgs 
doublets, and the magnitude of /i, but not its sign, is determined by the radiative electroweak breaking (REWSB) condition. In 
this paper we mainly focus on /i < as well as M 2 < 0, as explained earlier. Although not required, we will assume that the 
gauge coupling unification condition g 3 = g\ = g 2 holds at Mqut in 4-2-2. Such a scenario can arise, for example, from a 
higher dimensional 50(10) [14] or 5/7(8) [15] model after suitable compactification. 



III. PHENOMENOLOGICAL CONSTRAINTS AND SCANNING PROCEDURE 

We employ the ISAJET 7.80 package [16] to perform random scans over the parameter space listed in Eq.(3). In this package, 
the weak scale values of gauge and third generation Yukawa couplings are evolved to Mqut via the MSSM renormalization 
group equations (RGEs) in the DR regularization scheme. We do not strictly enforce the unification condition g 3 = g\ = g 2 
at Mqut, since a few percent deviation from unification can be assigned to unknown GUT-scale threshold corrections [17]. 
The deviation between g 1 = g 2 and g% at Mqut is no worse than 3.5%, and it is also possible to get perfect gauge coupling 
unification. Perfect gauge coupling unification is typically not possible with gaugino universality. With nonuniversality in the 
gaugino sector, one can adjust M 3 /M 2 appropriately to get exact gauge coupling unification. If neutrinos acquire mass via 
Type I seesaw, the impact of the neutrino Dirac Yukawa coupling in the RGEs of the SSB terms, gauge couplings and the third 
generation Yukawa couplings is significant only for relatively large values (~ 2 or so). In the 4-2-2 model we expect the largest 
Dirac Yukawa coupling to be comparable to the top Yukawa coupling (~ 0.6 at Mqut). Therefore, we do not include the Dirac 
neutrino Yukawa coupling in the RGEs. 

The various boundary conditions are imposed at Mqut and all the SSB parameters, along with the gauge and Yukawa 
couplings, are evolved back to the weak scale Mz. In the evaluation of Yukawa couplings the SUSY threshold corrections [19] 
are taken into account at the common scale Msusy = sf m i L m t R • The entire parameter set is iteratively run between Mz 
and Mqut using the full 2-loop RGEs until a stable solution is obtained. To better account for leading-log corrections, one- 
loop step-beta functions are adopted for gauge and Yukawa couplings, and the SSB parameters ra$ are extracted from RGEs at 
multiple scales mi = m^m^). The RGE-improved 1-loop effective potential is minimized at an optimized scale Msusy, which 
effectively accounts for the leading 2-loop corrections. Full 1-loop radiative corrections are incorporated for all sparticle masses. 

The requirement of REWSB [20] puts an important theoretical constraint on the parameter space. Another important constraint 
comes from limits on the cosmological abundance of stable charged particles [21]. This excludes regions in the parameter space 
where charged SUSY particles, such as f\ or ti, become the lightest supersymmetric particle (LSP). We accept only those 
solutions for which one of the neutralinos is the LSP and saturates the WMAP (Wilkinson Microwave Anisotropy Probe) dark 
matter relic abundance bound. 

We have performed random scans for the following parameter range: 
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< mo, rriH u , ^H d < 20 TeV 
-2TeV < M 2 < 
< M 3 < 2 TeV 
45 < tan/3 < 55 
-3 < A /m < 3 

M < 0, (4) 

with m £ = 173.1 GeV [22]. The results are not too sensitive to one or two sigma variation in the value of m t . We use 
^b(^z) = 2.83 GeV which is hard-coded into ISAJET. This choice of parameters was influenced by our previous experience 
with the 4-2-2 where we set /i, M 2 > 0. 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm as described in [23]. All of the collected data 
points satisfy the requirement of REWSB, with the neutralino in each case being the LSP. Furthermore, all of these points satisfy 
the constraint ^cdm^ 2 < 10. This is done so as to collect more points with a WMAP compatible value of cold dark matter 
(CDM) relic abundance. For the Metropolis-Hastings algorithm, we only use the value of ^cdm^ 2 to bias our search. Our 
purpose in using the Metropolis-Hastings algorithm is to be able to search around regions of acceptable ^cdm^ 2 more fully. 
After collecting the data, we impose the mass bounds on all the particles [24] and use the IsaTools package [25] to implement 
the following phenomenological constraints: 

m h (lightest Higgs mass) > 114.4 GeV [26] 
BR(B S < 5.8 x 1(T 8 [27] 

2.85 x 1(T 4 < BR(b sj) < 4.24 x 1(T 4 (2a) [28] 

0.15 < gff^^g^MSSM < / 3 x [28] 

VcDuh 2 = O.lllini? ( 5 ^) [29] 
3.4 x 10- 10 < A(g - 2) M /2 < 55.6 x lO" 10 (3a) [11] 

We apply the experimental constraints successively on the data that we acquire from ISAJET. 



IV. SIGN OF ii AND YUKAWA UNIFICATION 

To first appreciate the impact of the sign of fi on Yukawa coupling unification and to see why fi < is preferred over fi > 0, 
in Fig. 1 we show the evolution of the top, bottom and tau Yukawa couplings for a representative Yukawa coupling unification 
solution. We observe that Yukawa unification requires relatively large threshold corrections to y^. To quantify the magnitude of 
the threshold corrections, we scan the parameter space given in Eq.(4) and calculate the finite and logarithmic corrections to 5yi, 
where the index i refers to top, bottom and tau Yukawa couplings. 
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FIG. 1. Evolution of top (red), bottom (blue) and tau (green) Yukawa couplings for \i < 0. 
Following [3], we define the quantity R as, 
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max(y t ,y b ,y r ) 

= — —, r- P) 

min(2/ t ,2/ 6 ,2/ T ) 

Thus, is a useful indicator for Yukawa unification with R < 1.1, for instance, corresponding to Yukawa unification within 
10%. 

Fig. 2 shows a plot of 5yi/yi versus R for both \i > and /i < 0. Points in blue, green and red represent, respectively, dy^/y^, 
Sy T /y T and Sy t /y t . Note that we choose the sign of Sjji from the perspective of evolving yi from Mqut to Mz. Fig. 2 confirms 
the trend seen in Fig. 1 that y t and y T receive small threshold corrections compared to y^. Therefore, it is reasonable to focus on 
the threshold corrections to y b while studying Yukawa unification. 

The scale at which Yukawa coupling unification is to occur is set by gauge coupling unification and is Mqut- Let us first 
consider the case of ^(Mqut) ~ ^(^gut)- Because the threshold corrections to y t are very small, it is convenient to think 
of y T . The SUSY correction to the tau lepton mass Sm T is given by Sm T = v cos (35y r . In order to get the correct r mass (m r ), 
one has to get an appropriate 5y T . Because of the range of values of cos f3 for large tan /?, there is freedom to choose the value of 
Sy T . It may be possible to trade this freedom in favor of top-tau Yukawa unification ^(Mqut) ~ ^t(^gut)- One then needs 
the correct SUSY contribution to Sy^ in order to achieve Yukawa coupling unification ^(Mqut) ~ 2/6 (^gut) ~ y T (^gut)- 

In order to understand why sign of fi is very crucial for Yukawa unification condition lets first analyze the analytical expression 
of threshold corrections for the bottom Yukawa coupling. The dominant contribution to Sy^ comes from the gluino and chargino 



(a) fi > (b) fi < 




m (TeV) m (TeV) 



FIG. 3. Syb/yb as a function of mo- Points in black correspond to 10% or better Yukawa unification (R <= 1.1). 
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loops, and in our sign convention, is given by [19] 

x , .finite _ 93 MTOgton/? y 4 2 tan/3 

t 

where #3 is the strong gauge coupling constant, rrig is the gluino mass, m^ is the sbottom mass, m^ is the stop mass, and A t is 
the top trilinear coupling. One can see from Fig. 2 that in order to achieve Yukawa coupling unification R ~ 1, the threshold 
corrections to y\> have to be negative (in our sign convention for 5yi) and in a somewhat narrow interval (—0.5 < Sy^ jyh < — 1.5) 
considering the full range of possible values of Sy^. The logarithmic corrections to y^ are in fact positive. This leaves the finite 
corrections to provide for the correct Sy^ to compensate for the 'wrong' sign of the logarithmic corrections. If \i > 0, the gluino 
contribution is positive, and so the contribution from the chargino loop must cancel the contribution from the gluino loop and the 
logarithmic correction, as well as provide the correct (negative) contribution to Sy^. This can be achieved only for a large mo, 
as for large mo and for large A t , the gluino contribution scales as M^/wq while the chargino contribution scales as A t jm\. It 
also should be noted that the numerical factor for the gluino contribution is larger than than the corresponding factor for chargino 
contribution. Therefore, a sufficiently large value of A t and mo is needed. This large required value of A t is the reason behind 
the requirement of A /mo ~ —2.6 for \i > 0. 

The scenario with /i < is interesting because the gluino contribution to Sy^ has the correct sign to obtain the required 
b-quark mass. Thus, we should expect that with \i < 0, we can realize Yukawa unification for a wider range Aq values. With 
the threshold contribution to y T proportional to M 2 tan /3, and with M 2 as a free parameter in the 4-2-2 model, we also should 
expect Yukawa unification to occur over a broader range of tan fj values. 

To proceed further, in Fig. 3 we present Sy^/y^ versus mo for /i > and /i < 0, by performing a scan over the parameter 
space given in Eq.(4). Here Sy^ includes full one loop finite and logarithmic corrections, and the black points correspond to 10% 
or better Yukawa unification (R < 1.1). It is clear from Fig. 3 that Yukawa coupling unification with a relatively light mo ~ 400 
GeV can be realized for fi < 0. The /i > scenario typically requires a very large mo > 8 TeV. 

V. YUKAWA UNIFICATION AND SPARTICLE SPECTROSCOPY 

We now present the results of the scan over the parameter space listed in Eq.(4). In Fig. 4 we show the results in the R- mo, 
R - tan f3, R- Ao/mo and M3 - M 2 planes. The gray points are consistent with REWSB and Xi LSP. The blue points satisfy the 
WMAP bounds on Xi dark matter abundance, sparticle mass bounds, constraints from BR(B S — > and BR(b — >> 57). 

The green points belong to the subset of blue points that satisfy all constraints including (g — 2) M . In the M3 - M 2 plane, points 
in red represent the subset of green points that satisfies Yukawa coupling unification to within 10%. 

In the R - mo plane of Fig. 4 we see that with both /i < and M 2 < 0, we can realize Yukawa unification consistent with 
all constraints mentioned in Section III including the one from (g — 2) M . This is possible, as previously noted, because we can 
now implement Yukawa unification for relatively small mo(~ 400 GeV) values because \i < 0, and, in turn, (g — 2) M obtains 
the desired SUSY contribution which is proportional to /iM 2 . This is more than an order of magnitude improvement on the mo 
value required for Yukawa unification with /i > 0. We also see from the R - Ao/mo plane that, as explained earlier, Yukawa 
unification for an essentially arbitrary value of Aq is obtained for /i < 0. Our observation about relaxing the possible range 
of tan P that accommodates Yukawa unified models is explicitly shown in the R - tan f3 plane. This also suggests that we are 
perhaps somewhat conservative in limiting the range of tan /?. While it is beyond the scope of this paper, it would be nice to 
systematically search for a lower bound on tan (j that is consistent with third family Yukawa unification. 

The Ms - M 2 plane of Fig. 4 has some interesting features. The large gray regions appear because the relic density of 
neutralinos is too high. This may seem peculiar given that we have non-universal Higgs boundary conditions. However, it is 
readily explained by the fact that M 2 < 0. The bulk of the gray region has a neutralino that is too light (m^o < 30 GeV). While 
it appears possible to have a relic density of a relatively light neutralino (< 30 GeV) consistent with WMAP, BR(B S — >> 
and BR(b — > 87), so that some of this gray region may turn blue, there will always be regions where some of these constraints 
are not satisfied. 

In Fig. 5 we show the relic density channels consistent with Yukawa unification in the m~± - m^o, m^ - m^o, m? - m^o and 
mA - m^o planes. All of the points shown in this figure satisfy the requirements of REWSB, Xi LSP, particle mass bounds and 
constraints from BR(B S — ^ and BR(b — ^ 57). The light blue points satisfy, in addition the constraint from (g — 2) M . 

The green points form a subset of light blue points that satisfies Yukawa unification to within 10%. The orange points satisfy all 
the constraints mentioned in Section III, while the red points form a subset of orange points that have R < 1.1. This choice of 
color coding is influenced from displaying the sparticle spectrum with and without neutralino dark matter, while still focussing 
on all the other experimental constraints. The idea is to show the myriad of solutions that implement Yukawa unification and are 
consistent with all known experimental bounds except for the bound on relic dark matter density from WMAP. The appearance 
of a variety of Yukawa unified solutions with a very rich sparticle spectrum is a characteristic feature of fi < 0. 
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FIG. 4. Plots in the R - mo, R - tan (3, R - Aq/ttiq and M3 - M2 planes. Gray points are consistent with REWSB and xi LSP. Blue points 
satisfy the WMAP bounds on xi dark matter abundance particle mass bounds, constraints from BR(B S — ► BR(B U — >> tz/ t ) and 

— ► S7). Green points belong to the subset of blue points that satisfies all constraints including (g — 2) M . In the M 3 - M 2 plane, points 
in red represent the subset of green points that satisfies Yukawa coupling unification to within 10%. 



We can see in Fig. 5 that a variety of coannihilation and annihilation scenarios are compatible with Yukawa unification and 
neutralino dark matter. Included in the tua - rn^o plane is the line tua = 2m^o which indicates that the A funnel region is 
compatible with Yukawa unification. In the remaining planes in Fig. 5, we draw the unit slope line which indicates the presence 
of gluino, stau and wino coannihilation scenarios. From the m~± - m^o plane, it is easy to see the light Higgs (h) and Z 
resonance channels. Our results are focussed on relatively light neutralinos (m^o < 225 GeV). We expect, based on the 
discussion presented, that other coannihilation channels like the stop coannihilation scenario are also consistent with Yukawa 
unification but we did not find them because of lack of statistics. 



VI. YUKAWA UNIFICATION AND DARK MATTER DETECTION 

In light of the recent results by the CDMS-II [30] and Xenon 100 [31] experiments, it is important to see if Yukawa unification, 
within the framework presented in this paper, is testable from the perspective of direct and indirect detection experiments. The 
question of interest is whether /i ~ Mi is consistent with Yukawa unification, as this is the requirement to get a bino-higgsino 
admixture for the lightest neutralino which, in turn, enhances both the spin dependent and spin independent neutralino-nucleon 
scattering cross sections. In Fig. 6 we show the spin independent and spin dependent cross sections as a function of the neutralino 
mass. In the case of spin independent cross section, we also show the current bounds and expected reach of the CDMS and 
Xenon experiments. The color coding is the same as in Fig. 5. A small region of the parameter space consistent with Yukawa 
unification and the experimental constraints discussed in Section III (red points in the figure) is excluded as we can see, by the 
current CDMS and XENON bounds. This shows that the ongoing and planned direct detection experiments will play a vital role 
in testing Yukawa unified models. 
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FIG. 5. Plots in the ~± - m^o , rrig - m^o , m? - m^o and tua - m^o planes. All points satisfy the requirements of REWSB, \i LSP, particle 

mass bounds and constraints from BR(B S — »• BR(B U —> tv t ) and BR(b — > sj). Light blue points further satisfy the constraint 

on (g — 2)fj,. Green points form a subset of light blue points that satisfies Yukawa unification to within 10%. Orange points satisfy all the 
constraints mentioned in Section III while red points form a subset of orange points that have R < 1.1. 
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FIG. 6. Plots in the asi - m^o and ctsd- m^o planes. Color coding is the same as in Fig. 5. In the asi - m^o plane we show the current bounds 
(black lines) and future reaches (red lines) of the CDMS (solid lines) and Xenon (dotted lines) experiments. In the ctsd - m^o plane we show 
the current bounds from Super K (black line) and IceCube (dotted red line) and future reach of IceCuce DeepCore (red solid line). 



In the case of spin dependent cross section, we show in Fig. 6 the current bounds from the Super-K [32] and IceCube [33] 
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experiments and the projected reach of IceCube DeepCore. It should be noted that IceCube currently is sensitive only to relatively 
large neutralino masses and therefore does not constrain the parameter space that we have considered. Likewise, while Super-K 
is sensitive in this region, the bounds are not stringent enough to rule out anything. However, from Fig. 6 we see that the future 
IceCube DeepCore experiment will be able to constrain a significant region of the parameter space. 

It is interesting to comment on the correlation between spin independent and spin dependent cross sections. While both of 
these cross sections are enhanced by the presence of a larger higgsino component in the neutralino, the spin independent cross 
section falls off as \jm\. In Fig. 7 we plot the correlation of the spin independent cross section versus the spin dependent cross 
section for rn^o > 50 GeV. We impose this lower bound on the neutralino mass only because for 50 GeV < rn^o < 400 GeV, 
both the spin independent and spin dependent future reach and spin independent current bounds are nearly flat. This provides 
one with an opportunity to identify regions with simultaneously high spin independent and spin dependent cross section which 
may have the best hope for being tested in ongoing and future experiments. 

In Fig. 7 we also show lines corresponding to a si = 4 x 10 -8 , 1.5 x 10 _9 (pb) and osd = 1-5 x 10 _5 (pb). These lines serve 
as a guide to demonstrate the cross correlation of spin independent and spin dependent cross sections and the plausibility of either 
discovering a Yukawa unified neutralino dark matter solution or definitively ruling out certain regions of the allowed parameter 
space. It is to be noted that since ongoing and future experiments require large spin independent and spin dependent cross 
sections, neutralino coannihilation channels such as stau coannihilation, bino-wino coannihilation and gluino coannihilation 
may be difficult to test with direct as well as indirect detection experiments such as IceCube DeepCore. 

Let us remark on the low mass neutralinos that we have found in this model. Because of the relative negative sign between 
M2 and M3, it is possible in principle to have Mi ~ 0. This means that within this framework, the neutralino may be as light as 
we want. The neutralino mass nonetheless is bounded from below because of the relic density bounds on dark matter. The 4-2-2 
model as we have presented it has all the ingredients needed to lower the neutralino mass to the lowest possible value allowed by 
various constraints. We do not focus on finding the lightest neutralino in this model. Of the data that we collected, the lightest 
neutralino found that is consistent with Yukawa unification has mass ~ 43 GeV. If we do not insist on Yukawa unification, we 
can get a neutralino as light as ~ 32 GeV. If we are willing to give up neutralino dark matter, then Mi ~ consistent with 
Yukawa unification is possible as is evident in Fig. 5. This, of course, requires invoking some other dark matter candidate such 
as axino. 



T I I I 




(Pb) 



FIG. 7. Correlation of asi and ctsd- Color coding same as in Fig. 5. Also shown are lines corresponding to <tsi = 4 x 10 8 , 1.5 x 10 9 (pb) 
andasD = 1.5 x 10~ 5 (pb). 

Finally in Table I we present some benchmark points for the 4-2-2 Yukawa unified model with \i < 0. All of these points 
are consistent with neutralino dark matter and the constraints mentioned in Section III. Point 1 corresponds to a solution with 
perfect Yukawa unification (R = 1.0). Point 2 represents a solution with minimum neutralino mass (43 GeV) consistent with 
Yukawa unification. Point 3 has a significant bino-higgsino admixture and, therefore, has relatively large spin independent and 
spin dependent neutralino-proton scattering cross-sections. Point 4 depicts a solution with a stop mass of only 826 GeV. Finally 
in point 5 we show an example with bino-gluino coannihilation channel with the gluino as light as 259 GeV. This should be 
relatively easy to find at the LHC. 

VII. CONCLUSIONS 

We have shown that Yukawa coupling unification consistent with known experimental constraints is realized in a SUSY 
SU (4) c x SU (2)l x SU (2)r model. With fi < Yukawa coupling unification is achieved for mo > 400 GeV, as opposed to 
mo > 8 TeV for fi > 0, by taming the finite corrections to the b-quark mass. By considering M2 < and M3 > gauginos 
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and fx < 0, we can obtain the correct sign for the desired contribution to (g — 2) M . This enables us to simultaneously satisfy the 
requirements of t — b — r Yukawa unification, neutralino dark matter and (g — 2) M , as well as a variety of other known bounds. 
We have demonstrated the existence of a variety of coannihilation scenarios involving gluino, wino and stau, in addition to the 
light Higgs, Z and A resonance solutions. The Yukawa unified solutions may also have relativley large spin independent and 
spin dependent interaction cross sections with nucleons in the case of mixed bino-Higgsino dark matter. Finally, within the 4-2-2 
model, it is possible to obtain relatively low neutralino masses ~ 43 GeV (~ 30 GeV without Yukawa unification) consistent 
with neutralino dark matter. 





Point 1 


Point 2 


Point 3 


Point 4 


Point 5 


mo 


1UZ / 


loUU 


1Z1U 


yoK) 


1 /zu 


Mi 


-665 


-81 


-414 


-126 


-538 


Mo 


-1475 


-543 


-940 


-517 


-943 


M 3 


550 


611 


374 


460 


70 


tan j3 


49.1 


52.8 


50.6 


47.0 


47.6 


A /m 


0.26 


1.06 


-1.15 


-1.08 


-1.25 


rriHu 


743 


1919 


1231 


1090 


295 


rriHd 


1505 


2395 


1745 


1869 


1729 


m h 


114 


115 


114 


115 


115 


m H 


847 


573 


781 


1100 


1006 


m A 


841 


569 


776 


1090 


1000 


m H ± 


852 


581 


787 


1100 


1010 


^1,2 


280,341 


43,352 


168,242 


56,337 


233,782 




352,1236 


380,513 


246,795 


371,476 


1210,1216 




342,1225 


355,509 


239,786 


338,475 


782,1217 


m g 


1321 


1470 


955 


1110 


270 




1771,1489 


2170,2130 


1550,1410 


1400,1320 


1818,1697 




1053,1410 


1400,1440 


822,1040 


826,965 


1070,1248 


m d L , R 


1773,1512 


2180,2160 


1550,1440 


1400,1370 


1820,1730 


mr 

0l,2 


954,1399 


1350,1430 


774,1020 


724,906 


992,1245 




1391 


1810 


1340 


1000 


1807 




1211 


1420 


1100 


759 


1550 


m e L , R 


1393,1096 


1820,1820 


1340,1250 


1010,1040 


1809,1763 


171 T 1, 2 


500,1212 


885,1420 


641,1110 


462,765 


1170,1554 


CTSl(pb) 


4.02 x 10~ 8 


4.1 x 10~ 9 


4.1 x 10~ 8 


9.5 x lO -10 


1.1 x lO -10 


(pb) 


8.4 x 10~ 5 


7.5 x 10~ 6 


1.7 x 10- 4 


8.2 x 10~ 6 


2.9 x 10~ 8 


^CDMh 2 


0.08 


0.11 


0.09 


0.08 


0.11 


R 


1.01 


1.11 


1.09 


1.07 


1.08 


^i(Mgut) 


0.98 


0.98 


0.99 


0.98 


1.00 



TABLE I. Sparticle and Higgs masses (in GeV), with m t = 173.1 GeV. All of these benchmark points satisfy the various constraints 
mentioned in Section III and are compatible with Yukawa unification. Point 1 exhibits 'perfect' Yukawa unification, point 2 has the lightest 
neutralino, point 3 shows 'large' spin independent and spin dependent cross-sections, points 4 and 5 correspond to the lightest stop and gluino 
respectively. Point 5 also provides an example of bino-gluino coannihilation channel. 
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